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I. INTRODUCTION 

Trions (also called charged excitons) are bound states 
of a neutral exciton (electron-hole pair, X — e + h) with 
an additional carrier, either electron or valence hole, for 
a negative {X~ = 2e+h) or positive trion (X+ = 2h + e), 
respectively.^'^ Neutral and charged excitons occur nat- 
urally in photoluminescence (PL) experiments, in which 
creation or annihilation of e-h pairs accompanies inter- 
band absorption or emission of light »2. While excitons are 
expected in charge-neutral systems, formation of trions 
depends on the presence of free carriers. 

The trion binding energy A is defined as the effective 
attraction between an exciton and an additional carrier 
(A = E[X] + E[e] - E[X^], where E[. . .] is the ground 
energy of a given complex). An important material pa- 
rameter affecting A is the ratio of electron and hole ef- 
fective masses,-^ rj = fie/fih', for 77 = the X~ dynamics 
reduces to a familiar problem of the D~ centeri^^ Orig- 
inating from charge-dipole interaction, trion binding is 
usually much weaker than the exciton binding. However, 
its strong enhancement by spatial confinement and/or 
magnetic field B was predicted.- The pioneering exper- 
iment of Kheng et alwk in CdTe, followed by a series of 
measurements in GaAsSiSiiSiii and ZnSeiSiH confirmed 
that it is sufficient for the trion's detection. The Ar+ was 
also observed^*' and shown to be different from X~ 
(due to the <-> rj~^ asymmetry). 

In a typical experimental configuration, a quantum 
well containing a quasi-two-dimensional (2D) electron (or 
hole) gas is placed in the perpendicular field. The field 
should be sufficiently strong to induce magnetic quantiza- 
tion of single-particle states into macroscopically degen- 
erate Landau levels (LLs), with the inter-LL (cyclotron) 
spacing TilOc which is at least comparable to the effective 
excitonic Rydberg {Ry « 5.5 meV in GaAs). The well 
width w should not be much greater than the effective 



excitonic Bohr radius [qb ~ 10 nm in GaAs) and the 
magnetic length (A = \fhcfeB w 8.1 nm at i? = 10 T). 

The quantum-mechanical problem of a 2D exciton in 
a high magnetic field goes back several decades liSiiLiS 
The continuous (owing to the neutral charge) energy dis- 
persion in some idealized situations is known analyti- 
cally, and the deviations due to inter-LL mixing or fi- 
nite well width have been studied experimentally and 
numerically»i2i The simple 9 = optical selection rule 
(g being the wavevector of an annihilated e-h pair) 
results from the small momentum carried by the ab- 
sorbed/emitted photon. In experiment, the excitonic re- 
combination is usually observed in PL of quantum wells 
containing no free carriers. 

2D trions at high B have also been extensively stud- 
ied in the past.^" Different theoretical/computational ap- 
proaches and the key results have been discussed in an ex- 
haustive review by Peeters, Riva, and Vargai2i Depend- 
ing on the parameters (e.g., w and i3), the trion energy 
spectrum contains one or more bound states, which can 
be conveniently labeled by total spin S of the pair of 
electrons and total angular momentum M . At small -B, 
the only bound X^ state is the spin-singlet with S" = 
and M — 0, equivalent to a 2D Hydrogen iom^iSiSiifliii 
In the (unrealistic) limit of very high B and vanishing w, 
the singlet unbinds, and it is replaced by the spin-triplet 
with 5=1 and M = — li^^i^^ A triplet trion was iden- 
tified by Shields et alm& (although now it is not clear if 
it was the M = —1 triplet ground state). The singlet- 
triplet crossover was predicted^i in rather high fields (e.g. 
B « 30 T in narrow symmetric GaAs wells) . Despite ini- 
tial difficulties^^ it was also eventually confirmed in sev- 
eral experimentsi2Si2L22i2S Additional bound states were 
predicted^ at intermediate B, but (for realistic parame- 
ters) they are always less strongly bound than the above 
two. These states were also confirmed by both indepen- 
dent calculations'^^ and by experimentsi^2i22iS 
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The pair of optical selection rules results from in- 
variance under (magnetic) translationa^^*^ and from 
the particle-hole symmetry between conduction electrons 
and valence holes i^^i^SiSi Both these "geometric" and 
"hidden" symmetries are at least weakly broken in re- 
alistic conditions. Nevertheless, recombination of trions 
with M ^ requires a symmetry-breaking collision (with 
an impurity or another carrier), and therefore, it became 
customary to label different trions as "bright" and "dark" 
(meaning having M = and M ^ 0, respectively). In 
most recent experiments up to three trion states 
are observed: "bright" singlet X~ , "dark" triplet 
and a weakly bound "bright" triplet with S = \ 
and M = 0. Vanishing oscillator strength, of X^^ 
confirmed directly in optical absorption,'''* makes it more 
difficult to observe (even in PL) than the other trions. In 
fact, this state had not been conclusively identified un- 
til the work of Yusa et al.r^ motivated by earlier high- 
field experiments (especially of Hayne et al.^ showing 
apparent discrepancy with prediction— of singlet-triplet 
crossing, but also several others2& showing inconsistent 
features in trion spectra) and a numerical prediction"''^ 
of an excited triplet state, X~^. Most recently, ad- 
ditional weakly bound states were reportedSS in CdTe 
{Ry ~ 10 meV, about twice larger than in GaAs). 

Involving only three particles, trion's quantum dynam- 
ics might appear to be relatively simple, both concep- 
tually and computationally. Addition of high magnetic 
field and confinement does not add much complexity, and 
indeed, quantum numbers and symmetries of all bound 
trions (in wide range of realistic conditions) are by now 
understood. In short, the dynamics depends on compe- 
tition of several energy scales, with the following char- 
acteristic values for our example of a 15 nm symmetric 
GaAs weU at _B = 10 T: cyclotron gap {TiuJc ~ 18 meV 
for electrons and ~ 4 meV for heavy holes). Coulomb 
energy (e^/A ~ 14 meV), and a small Zeeman gap. The 
well width w comes in two places, defining excitation 
gaps to higher subbands (gap to the second subband is 
~ 50 mcV for electrons and ~ 10 meV for the holes) and 
affecting in-plane Coulomb matrix elements (this effect is 
parametrized by w/\ ~ 2). Also, even weak asymmetry 
between electron and hole subband wavefunctions Xs{z) 
can considerably affect the trion binding, since it leads to 
different magnitudes of e-e and e-h interactions, which 
no longer cancel in A (this effect is essential in asymmet- 
ric wells, not considered here, where the asymmetry de- 
pends on carrier concentration and/or additional gates). 

However, the above optimistic (and popular) view 
hides the fact that good understanding of the role of tri- 
ons in PL experiments must include variety of complica- 
tions due to coupling to the environment (e.g., complex 
single-particle energy band structure with nonparabolic 
and anisotropic dispersions, interaction with free car- 
riers, lattice defects, and phonons, spin-orbit effects, 
etc.). Some quantities (e.g., the binding energies, es- 
pecially of the triplet states) can be calculated rather 
accuratelyi^S*^ However, others (such as the critical val- 



ues of B and w for the singlet-triplet transitions) depend 
so sensitively on the (often unknown) parameters that 
their quantitative modeling turns out somewhat point- 
less. Another unsolved (quantitatively) and important 
problem is the kinetics of trionS)22ii£ involving their bind- 
ing/unbinding [X^ X + e~) and, at high _B, or- 
bital/spin relaxation [X^ ^ X^^ ^ ^tb)- 

In this article, we analyze the effect of trions on the op- 
tical absorption spectrum of 2D electrons in a magnetic 
field. Thus, in addition to the earlier studied bound trion 
statespS*^ the entire low-energy 2e + h spectrum is con- 
sidered. In a somewhat related work, asymmetry of trion 
absorption peaks at small B was discussed by Stebe et 
alMi The inclusion of only two electrons and one hole in 
the model restricts it to the "dilute" regime, defined by a 
small value of the filling factor, ^ ^ 1 (^ is defined as the 
number of electrons N divided by the LL degeneracy g; 
alternatively, i' — 2TTgX^ where g is the 2D electron con- 
centration). It is remarkable that in some systems the 
immersed trion is only weakly (perturbatively) affected 
by the surrounding electrons. In narrow wells and in high 
magnetic fields, this occurs at < ^ due to "Laughlin 
correlations" ^? between electrons and trions?^ prevent- 
ing strong e-X" collisions (indeed, in wider wells trions 
seem to be strongly coupled to the electrons and cannot 
be regarded as simple three-body quasiparticles'*^) . 

The density of states (DOS) and absorption spectra are 
calculated numerically by exact diagonalization in Hal- 
dane's spherical geometry^^i^ The figures were drawn 
for a particular choice of a symmetric w = 15 nm GaAs 
quantum well at i? = 10 T. Since absorption into bound 
trion states was studied earlier, we concentrate on the 
effects of the interaction of the exciton with a (single) 
unbound electron. The main conclusions are redistribu- 
tion of the density of states (by filling the gaps between 
the LLs) and the emergence of additional, trion peaks in 
the absorption spectrum (also in the excited LLs). 

The calculation is supplemented with the results 
of experimental polarization-resolved photoluminescence 
(PL) and photoluminescence-excitation (PLE) studies of 
a 2D hole and electron gases in a symmetric 15 nm GaAs 
well. Presented spectra reveal absorption into a pair of 
bright trions in the lowest LL and emission from these 
two states, the exciton, and the dark triplet trion. 



II. MODEL 

The e + h and 2e -\- h energy spectra are calculated 
by exact diagonalization of the hamiltonian matrix on 
a Haldane sphere, convenient for modeling an infinite 
plane with 2D translational invariance. In this geom- 
etry, the magnetic field normal to the spherical sur- 
face of radius R is produced by a Dirac monopole od 
strength 2Q = AnR? /(fiQ (^o = hc/e is the magnetic flux 
quantum). Through the following relation, R^ = QX^, 
monopole strength determines surface curvature (in the 
magnetic units). 
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The single-particle states are called monopole 
harmonics 4^ They are the eigenstates of angular momen- 
tum I and its projection m on the z-axis. The lowest shell, 
corresponding to the lowest LL, has I = Q and finite de- 
generacy g = 2Q + 1. Higher shells, corresponding to the 
excited LLs labeled by index n, have I = Q + n. 

The energy of the nth shell is e„ — Tiwdn + \) + 
'h'^n{n + l)/2fiR^. It contains the cyclotron gap (with 
Uc — eB / ^c, where /i is the effective mass) and an ad- 
ditional, curvature-dependent term. In order to model 
a real structure, we take advantage of the LL structure 
of monopole harmonics, but replace £„ with the correct 
energies of electron and hole LLs. E.g., at B = 10 T, we 
use hujc = 17.8 meV for electrons and 3.7 mcV for the 
heavy hole (the latter value taken after Cole et alM- for 
w = 15 nm). 

The second-quantization hamiltonian reads 



^c|'c]cfcQWyfei- (1) 
ijkl 



An additional term V cJctM,,- will also be included to 
describe interaction with a positive or negative point 
charge. Here, c| and Ci are operators creating and anni- 
hilating a conduction electron or a valence hole, in the 
state labeled by a composite index i containing all rele- 
vant single-particle quantum numbers (band /?, subband 
s, LL index n, angular momentum m, and spin a). 

The Coulomb interaction matrix elements v and u in 
the basis of monopole harmonics are known analytically 
for all particles confined to a 2D surface of the sphere. 
However, to account for a finite width of the real quantum 
well, we have integrated all matrix elements numerically, 
in 3D, using form- factors appropriate for the actual elec- 
tron and hole subband wavefunctions Xs{z)- Mixing with 
higher quantum well subbands (excited states in the z- 
direction, labeled by s > 0) is not very strong in a narrow 
and symmetric well due to high quantization energy (for 
w = 15 nm and Alo.35Gao.65As barriers, it is 49.1 meV 
and 11.5 meV to the first excited electron and hole sub- 
bands, respectively; calculation after Ref. fisf) and the 
parity conservation. Nevertheless, it is not quite negligi- 
ble for the bound states (states with strong interactions). 

Without impurities, the e + h and 2e + h eigenstates 
of H are labeled by total angular momentum L and its 
projection L^. When converting these quantities to the 
planar geometry, neutral and charged states must be 
treated differently: L of an e-h pair represents wavevec- 
tor q — L/R, and for a, 2e + h state it corresponds to 
M = L — Q. The 2e + h eigenstates are also labeled by 
spin S of the pair of electrons and its projection Sz- The 
calculation need only be performed in the = Sz = 
subspace and the appropriate Zeeman shift can be added 
at the end to the energies of each triplet {S = 1) state. 

With an impurity placed at a north pole of the sphere, 
Lz is still conserved, but L is not. Only the Sz = sub- 
space need be considered, but a separate diagonalization 
must be performed for each Lz- 
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FIG. 1: (color online) Density of states (DOS) of e + h and 
2e + h systems in a symmetric GaAs quantum well of width 
1(7 = 15 nm at magnetic field _B = 10 T, calculated on Haldane 
sphere with a large magnetic monopole strength 2Q = 50 in- 
cluding only two lowest electron and hole LLs {n < 1). Con- 
tinuous curves were obtained from finite-size discrete energy 
spectra by gaussian broadening. For 2e+h, DOS is divided by 
g = 2Q + 1 (LL degeneracy of the second electron), and singlet 
and triplet configurations (spin 5 = and 1) are drawn sepa- 
rately. For e + h, LL indices "n^nh" mark the strongest peaks. 
Inset: magnified lowest-energy sector (vertical lines mark po- 
sition of discrete bound trions and the excitonic ground state. 



The 2e + h diagonalization was carried out in 
configuration-interaction basis, \i,j;k) = cj^ctc^ | vac). 
Here indices i and j denote the occupied electron states, 
k describes the hole, and |vac) is the vacuum state. Sim- 
ilarly, the basis for the e + h calculation was \i;k). As 
mentioned earlier, only the spin-unpolarized states with 
Lz ^ rrii + nij — nik = are included in the basis. To 
find eigenstates corresponding to given (L, S) we used a 
modified Lanczos algorithm, with additional projection 
of Lanczos vectors at each iteration. 



III. RESULTS AND DISCUSSION 

A. Density of States 

We begin by the calculation of the density of states 
DOS = dT/dE (the number of states F per unit of en- 
ergy E) . In Fig. n we compare the results for e + h and 
2e + h obtained for a rather large 2Q = 50 and including 
only the lowest two electron and hole LLs {n < 1). The 
discrete energy spectra obtained from finite-size calcula- 
tion were converted to the continuous lines shown in the 
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figure by broadening witli Gaussians, 



Dos(£;) ^^^^(i;-^;,;), 



(2) 



where the summation goes over ah energy levels Ei and 
gg(x) = 5-^-1/2 exp(-x^/52). Thin and thick lines 
correspond to two different broadening widths S = 0.2 
and 0.5 meV. Blue solid, red dashed, and black dotted 
lines were used to distinguish 2e + h {S = and 1 plot- 
ted separately; only 5'^ = is shown for 5 = 1) from 
e + h. To compare DOS for e + h and 2e + h, the latter 
curve was rescaled from the value defined by Eq. ||2Jl by 
= {2Q + (i.e., divided by the number of states 
available to the second electron). Energy E is measured 
from the noninteracting configurations in the lowest LL. 

For a noninteracting e—h pair, DOS consists of dis- 
crete peaks labeled by the LL indices for both particles, 
"ngn/i". The black dotted curve demonstrate the effect 
of e-h interaction. Within each pair of LLs, magneto- 
excitonic dispersion becomes flat at long wavevectors q, 
corresponding to the vanishing e-h attraction. There- 
fore, although smeared toward lower energies, strong 
^^neUh" peaks persist. In our calculation, restricted area 
of the sphere prohibits the pair to separate in the q oo 
limit (the range of q = L/R is restricted by L < 2Q), 
and the "Uerih^ peaks have finite height and are dis- 
placed to the left by the remnant attraction (^ e^/2R — 
(2yQ)-i eVA). This redshift (- 1.4 meV for 2Q = 50) 
is a finite-size artifact. The low-energy tail extends from 
each "rien/i" peak over the range of the excitonic bind- 
ing energy within the corresponding pair of LLs. Since 
e^/A is larger than hujc of the holes, these tails essen- 
tially close the gaps between the neighboring hole LLs. 
Thin lines (shown in the magnified low-energy sector in 
the inset) reveal artificial size quantization on a sphere 
(i = 0, 1, 2, . . . ). On the other hand, an interesting real 
feature is the maximum at the beginning of the "10" tail 
{E w 10 meV). It is due to the fact that excitonic disper- 
sion for rig 7^ Uh is nonmonotonic at small wavevectors q 
(the ground state occurs at g > 0, leading to dE/dq = 
and dT/dq cx g > 0, i.e., to a singularity in dT/dE). 

The curves for 2e -\- h are also dominated by the non- 
interacting peaks ^^Uen'^rih" corresponding to three un- 
bound particles confined to different combinations of LLs. 
Low-energy tails describe the exciton with an additional 
unbound electron, (note oscillations due to excitonic size 
quantization). Bound trion states appear as discrete 
peaks below the continuous tails, well visible only in the 
inset, additionally marked with red and blue bars. All 
three trions: X~ , X~^, and appear bound (excitonic 
ground state is marked by a black bar for comparison). 
However, weak (< 1 meV) binding of , virtually equal 
to X~^, is an artifact of the unrealistic n < 1 restriction. 

The emergence of bound trion states below the exci- 
ton's continuum is one noticeable difference between the 
e + h and 2e + h DOS. Another significant effect of the 
exciton-electron interaction is further (compared to one 
due to excitonic e-h interaction) smearing of the LLs, 
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FIG. 2: (color online) Binding energies A of four different 
negative (a) and positive (b) trions (the pair of labels 
in parentheses are two-electron spin S and total angular mo- 
mentum M) in a symmetric GaAs quantum well of width 
w — 15 nm as a function of magnetic field B, calculated on 
Haldane sphere with a magnetic monopole strength 2Q — 20 
including five lowest LLs (n < 4) and either one (s = 0) 
or two lowest subbands (s < 1) for each electron and hole. 
Continuous curves were obtained by interpolation from exact- 
diagonalization at the values of B marked with dots. 



i.e., transfer of the density of states away from LLs and 
filling the gaps between them. Also, DOS within the 
inter-LL regions shows features related to the interac- 
tions in a nontrivial manner. For example, spin depen- 
dence of structures at i? « 5 meV and E' « 25 meV re- 
veals their connection with the (obviously, spin-sensitive) 
exciton-electron interaction. Note that the e + h peak at 
E 10 meV, identified earlier with an inter-LL exciton, 
persists in the 2e + h curves regardless of spin configura- 
tion (to confirm its one-electron nature). 

The n < 1 restriction to only two lowest LLs was help- 
ful in demonstrating LL smearing and emergence of ad- 
ditional peaks between LLs due to e~h and X-e inter- 
actions used in Figs. 1 and 2. However, it gives incor- 
rect exciton and trion energies and, more importantly, ig- 
nores the contribution to DOS coming from the neglected 
higher hole LLs (recall that hole's hcuc is only 3.7 meV 
at i? = 10 T, much smaller than electron's 17.8 meV). 

More accurate trion binding energies (for both X^ and 
X+) are plotted in Fig. |21as a function of B. These val- 
ues were obtained by including five LLs {n < 4) and up 
to two subbands (s < 1) for each electron and hole. The 
lowest-subband (s = 0) calculation is similar to Ref. Isol 
but it used more accurate Coulomb matrix elements, cal- 
culated for the actual Xo{z)- For X^, an additional. 
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FIG. 3: (color online) The same as in Fig.[Tlbut for 2Q = 20, 
n < 4, and without dividing 2e + h DOS hj g = 2Q + \. 



weakly bound "dark singlet" state occurs only at quite 
high fields. The curves for s < 1 demonstrate that in a 
narrow symmetric well the subband mixing is most ef- 
ficient for the M — Q singlet, while the lowest-subband 
approximation is accurate for both triplets. 

In our calculations we assumed the Al fraction of 
X — 0.35 in the Ala;Gai_a;As barriers. This restriction 
is justified by the fact that the binding energies of 
are almost insensitive to the increase of Al fraction all 
the way up to a; = 1. The only exception is data for X+ 
marked with open dots in Fig. |2Ib), obtained for pure 
AlAs barriers. Evidently, the A"+ binding energies de- 
pend more strongly on the barrier height, which must be 
taken into account in a realistic model. 

More accurate DOS is shown in Fig. O including five 
LLs in < 4) but still only the lowest subband. Since 
dimension of the Hilbert space quickly grows with ^max 
and 2Q, we were forced to use a smaller 2Q — 20 in 
this case (for 2e + h this yields dimension 58,875 for the 
Lz = Sz = subspace, whose all eigenenergies must 
be calculated to plot DOS). Especially at higher ener- 
gies, the curves would become very complicated due to 
an increasing number of overlapping peaks correspond- 
ing to different combinations of rie and Uh - if not only 
5, but all LLs were included for the hole. By count- 
ing electron cyclotron gaps from the lowest LL peak at 
E Ri ~e^/2R, energies corresponding to consecutive n^s 
have been marked with gray vertical lines. Comparison 
of the curves for e + h (here plotted without rescaling by 
g^^) and 2e + h shows that the effect of smearing the 
LLs and filling the gaps between them due to X-e in- 
teraction is only enhanced when more LLs are included. 
In the inset, the binding energies of all trion peaks are 



FIG. 4: (color online) Density of states (DOS) of e -f /i in 
a symmetric GaAs quantum well of width w — 15 nm at 
magnetic field B = 10 T, in the presence of a positive (a) 
and negative (b) impurity at different distances d from the 
center of the well, calculated on Haldane sphere with magnetic 
monopole strength 2Q — 20 including five LLs (n < 4) and 
three subbands (s < 2) for electron and hole. Continuous 
curves were obtained by gaussian broadening. Marked bound 
states: D° ^ D+ + e, D+X = D+ + X, A° ^ A- + h, and 
A~ X — A" + X; e* and h* denote electron and hole in the 
excited (n = 1) LL. 



already well converged for n < 4 (note however that the 
lowest subband approximation considerably affects espe- 
cially the singlet state); the exact values are A = 1.72, 
0.93, and 0.24 meV for X~, X~^, and X^^, respectively. 

An ionized impurity can have a similar effect on the 
e+h DOS to that of an additional electron. This impurity 
can be either a positive or negative point charge (ionized 
donor Z?+ or ionized acceptor A~) placed at a distance d 
away from the center of the quantum well. The effects of 
and A~ are not equivalent due to the e-h asymmetry. 
Two frames of Fig.^show how the e + h DOS changes in 
the presence of D+ or A~ placed at d/w — ^ and (edge 
and center of the well). In plotted energy range (up to the 
first excitonic maximum "00"), the curves for d/w = 1 
are almost identical to those without an impurity). One 
or more bound D~^X or A~X states (analogous to trions) 
emerge below the excitonic tail, at the position sensitive 
to d. Being localized (and nondegenerate), they do not 
contribute to the macroscopic DOS. 

When the impurity approaches the well, a strong peak 
detaches from "00" and moves to the left through the 
excitonic tail. It corresponds to the (macroscopically de- 
generate) + h or A^ + e configuration with the un- 
bound e~h pair. In our example, it passes the tail's edge 
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{X ground state) when the impurity is aheady inside the 
well, at d « 3 nm for D"*" and 5.5 nm for A~ . For d — 0, 
especially the A" is bound much more strongly than X. 
Thus, the strongest effect of an impurity is that, depend- 
ing on d, it can bring down macroscopic DOS correspond- 
ing to unbound e-h below the free excitonic ground state. 
Certainly, the localized D^X and A^X states involving 
the e—h binding still remain the absolute ground states. 
Here, the qualitative difference caused by the impurity 
is that DOS raises abruptly from essentially zero at the 
nondegenerate bound state to the continuum of degener- 
ate unbound states (instead of through an excitonic tail). 
Especially the marginally bound A~'X is hardly distin- 
guished from the continuum. 

For d = we have also marked the + h* and A'^ + e* 
peaks corresponding to the unbound hole or electron in a 
higher LL, and thus originating from the "01" and "10" 
peaks without an impurity. These peaks are separated 
from -I- h and A° -I- e by a cyclotron gap; thus, an 
impurity mixes the LL peak structure with an excitonic 
tail due to e-h interaction. 



B. Oscillator Strength 

Let us now turn to the calculation of oscillator strength 
n for the vac e + ft, and e <-> 2e -I- ft optical transitions 
("—!■" for absorption; "<— " for emission). For a pair of 
initial and final states, e.g., — e and (jj — 2e + h, it is 
calculated from Fermi's golden rule. 



E 

k 



144 



(3) 



Here, k = n,m, c] and k = n, — m, — a], and the 
summation runs over all electron states k in the conduc- 
tion band and all corresponding hole states k in the va- 
lence band. Note that according to convention of Eq. 

creates electrons or holes, depending on band index /3, 
leading to the reversed sign of m and a in k. The oscil- 
lator strength for the recombination of initial (f) = e + h 
or 2e + h states, in the latter case summed over all final 
■0 = 6 states, can be expressed as a function of the initial 
energy E = E^, 



(4) 



This is the e + h oi 2e + h "optical density of states" 
(ODOS). Alternatively, oscillator strength (if necessary, 
weighted by the occupation function for the initial 
states ^ — e) can be expressed as a function of the photon 
energy £ ^ E^ - E^, 



^{£) — il^^Q.^ S{Ecj, — Eji, — £). 



(5) 



This is the cj) absorption spectrum (equivalent to 0's 
ODOS for ip — vac, but not for ijj = e). 



to 
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empty dashed bars - no interaction 
empty solid bars - isolated LLs 
filled bars - mixed LLs 

green dotted curve - DOS 



cyclotron gap 
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FIG. 5: (color online) Optical density of states (ODOS) of 
e + h, equivalent to absorption spectrum of an exciton, in a 
symmetric GaAs quantum well of width w — 15 nm at mag- 
netic field _B = 10 T, calculated on Haldane sphere with a 
magnetic monopole strength 2Q — 20 including five lowest 
electron and hole LLs (n < 4). Results without e-h interac- 
tion, with interaction within isolated LLs (peaks marked with 
LL indices "nenh"), and with all interaction effects (includ- 
ing both intra- and inter-LL scattering) are shown. The e + h 
DOS from Fig.j^is drawn for reference. 



In Fig. O we plot ODOS of e + ft. The parameters 
{w — 15 nm and _B = 10 T) and Hilbert space (n < 4 
and 2Q = 20) are the same as in Fig. |3| If the Coulomb 
energy e^/A were much smaller than the cyclotron gaps, 
then the only optically active states would be the q — 
excitons with electron and hole confined to the same LLs 
{ue — rih = n). As shown with empty bars, these exci- 
tonic peaks "nn" in ri(f ) would have equal height and 
occur at E{n) = nh{ujc,e +^c,h) — Ax(n), where Ax('t-) 
is the exciton binding energy in the nth LL. These ener- 
gies coincide with the low-energy edges of the excitonic 
tails of the "nn" peaks in DOS (replotted with a thin 
dotted line), where degeneracy is absent and DOS van- 
ishes. Since the exciton binding Ax{n) decreases as a 
function of n, the separation between consecutive peaks 
also decreases and it is always larger than the cyclotron 
gap between the corresponding "nn" maxima in DOS. 
For w = 15 and B = 10 T, distances between peaks 
"00", "11", and "22" (calculated excluding LL mixing 
by setting Vijki = unless rii = rij = rik = ni = n) are 
Ell - Eqo = 24.7 meV and £^22 - En = 22.8 meV, both 
considerably larger than ft(wc,e + Wc,/i) = 21.5 meV. 

In reality, inter-LL scattering mixes g = states with 
different n. As shown with filled bars, this causes shifting 
of the peaks (and a further, small increase of the spac- 
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FIG. 6: (color online) Optical density of states (ODOS) of 
2e+h, in a symmetric GaAs quantum well of width w = 15 nm 
at magnetic field _B = 10 T, calculated on Haldane sphere 
with a large magnetic monopole strength 2Q — 50 including 
only two lowest electron and hole LLs (n < 1). Singlet and 
triplet initial 2e + h spin configurations (S = and 1), and 
ground and excited final e states {n = and 1) are drawn 
with different lines. Strongest peaks are marked by LL indices 
of the initial and final states, "n'/rienji" with rie = rih- Weak 
peaks associated with bound trions {X~) are also indicated. 
The 2e + h DOS from Fig. is drawn for reference. 



ing between neighboring peaks) and transfer of oscillator 
strength from higher to lower energy. Nevertheless, the 
effect is perturbative and consecutive peaks can still be 
labeled by n. For w = 15 and i? = 10 T, the first two 
gaps in f2 (calculated including all matrix elements Vijki) 
increase to En — Eqq = 26.4 meV and E22 — En = 
23.0 meV. The relative magnitudes of the lowest three 
peaks are fin/fioo = 0.55 and 1^22/^^00 = 0.46. 

In Fig. El we plotted ODOS of 2e + /i calculated us- 
ing the same Hilbert space (n < 1 and 2Q — 50) as in 
Fig. n Also for this larger system there is no apparent 
correlation between the features in DOS and ODOS. This 
follows immediately from the q = selection rule, equiva- 
lent to the requirement of nondegenerate relative motion 
of the recombining e-h pair, preventing high DOS of a 
2e + h state involving such a pair. 

In analogy to e + h, the main peaks can be labeled by 
"n'/ nn" , and they correspond to a g = exciton cre- 
ated/annihilated on the nth LL, in the presence of the 
second, ?/;-electron on the n'th LL. In those main peaks 
"n'/nn", the second electron is not bound to the cre- 
ated/annihilated exciton; they describe excitonic optical 
processes, weakly affected by the exciton-electron scat- 
tering. Optically active bound trions (with electrons and 
holes in different LLs) appear in form of weaker peaks dis- 



o 
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FIG. 7: (color online) The same as in Fig.Elbut for 2Q = 20 
and n < 4. The strongest peaks are identified as excitonic 
^'n' /nri" transitions or trion recombination {X~). 



placed from ^^nl jnvb'' by the binding energy. In the lowest 
LL, the only well resolved trion is the singlet (dark triplet 
by definition has = and bright triplet is too weakly 
bound to be distinguished from exciton in this energy 
scale). In the triplet trion at E « 2.5 meV one of the 
electrons is in the n = 1 LL. Although unstable against 
inter-LL relaxation, this state can form by capturing a 
"00" photo-exciton by a thermally excited electron. 

Note that "shake-up" transitions^ (n'/nn ^ n" with 
n' n", i.e., combination of "nn" recombination with 
n' ^ n" cyclotron excitation of the second electron) are 
forbidden^ for an isolated trion due to invariance under 
2D (magnetic) translations. This selection rule does not 
preclude replicas of an exciton and an unbound second 
electron, but these transitions have negligible intensity 
for V ^ g^^ <C 1 and cannot be identified in Fig. El 

Fig. □ is similar to Fig. El but it shows 2e + h ODOS 
calculated with 5 LLs taken into account (n < 4 and 
2Q = 20). The picture becomes fairly complicated, but 
the idea is the same. Dominant peaks correspond to ex- 
citonic transitions and can be labeled by "n'/nn" (in 
the calculated spectrum, the assignment is straightfor- 
ward due to angular momentum conservation in the op- 
tical transition, here causing the L = Q + n' selection 
rule) . A great number of weaker transitions that involve 
exciton-electron interaction emerge around the excitonic 
peaks. In particular, bright trions appear below "0/00" 
and "1/00" (singlet and triplet, respectively). 

Fig. jSl presents the 2e + h absorption spectra ^{£), 
calculated assuming that an electron in the initial state 
is cither in the lowest LL (ground state) or in a higher. 
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FIG. 8: (color online) 2e + h absorption spectra in a sym- 
metric GaAs quantum well of width w = 15 ran at magnetic 
field _B = 10 T, calculated on Haldane sphere with a magnetic 
monopole strength 2Q = 20 including five lowest electron and 
hole LLs (n < 4). Continuous curves were obtained by gaus- 
sian broadening. Singlet and triplet 2e + h spin configurations 
(5 = and 1) are drawn separately. Frames (a) and (b) corre- 
spond to the initial-state electron in the lowest or first excited 
LL (n' = or 1). The strongest peaks marked as "nn" corre- 
spond to excitonic absorption in the nth LL {e + h absorption 
spectrum from Fig. |K|is shown with gray bars for reference). 
Insets show magnified regions around the lowest two peaks 
"00" and "11". Trion peaks {X~) are identified. 
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FIG. 9: (color online) Polarized photoluminescence-excitation 
(PLE) spectrum measured in a symmetric w = Vb nm GaAs 
quantum well with hole concentration q — 1.5 ■ 10^^ cm~'^, in 
magnetic field _B = 14 T, at temperature T = 4.2 K. 



be seen most clearly in the insets, in which the vicinities 
of peaks "00" and "11" have been magnified. In reality, 
their intensity relative to the excitonic peaks will depend 
on the filling factor and can be much higher than in our 
model (which represents a very dilute system with only 
one free electron per g = 2Q + 1 states of the lowest LL). 
Remarkably, the spin of strong trion-related features cor- 
relates with the parity of n' — n: singlet {S = 0) peaks 
appear for n' ~ below "00" and for ti' = 1 below "11", 
while triplets (S = 1) occur for the opposite combina- 
tions of n' = with "11" and n' = 1 with "00". 



IV. EXPERIMENT 



n — 1 LL (due to thermal excitation) . To observe super- 
position of many closely spaced small peaks, each discrete 
e 2e + h transition was broadened with a gaussian of 
width S — 0.5 meV (main frames) or 0.2 meV (insets). 
As a reference, the vac —^e + h excitonic spectrum is 
shown with gray bars. 

The main result is that when the (great number of) 
main ODOS peaks "n'/nn" are shifted by — n'Tiujce 
to convert ^1{E) into ^{£), they all fall exactly onto ab- 
sorption spectrum of a bare exciton. Presence of an ad- 
ditional electron does not cause shifting or splitting of 
these main absorption lines, and it has insignificant ef- 
fect on their relative intensities. This demonstrates that, 
somewhat surprisingly, bare excitonic absorption is un- 
affected by a dilute electron gas (neither by renormaliza- 
tion of energy nor by transfer of intensity between LLs) . 
This result is obtained for a realistic quantum well, with 
significant electron-hole asymmetry and LL mixing. 

The effect of free electrons is emergence of additional 
(compared to the bare exciton) features in absorption 
spectrum, the strongest of them associated with the for- 
mation of trions. In Fig. |S1 trion absorption peaks can 



To verify some of presented calculations we per- 
formed polarization-resolved photoluminescence (PL) 
and photoluminescence-excitation (PLE) experiments on 
a symmetric w = 15 nm GaAs/AlAs quantum well, 
containing a valence holes gas with concentration g = 
1.5T0^^ cm~^ and subject to a strong magnetic field. The 
PLE experiment consisted of measuring polarization- 
resolved PL (emission intensity / as a function of emis- 
sion energy fout) for a series of excitation energies E-m- 
The spectra presented in Fig. El were recorded at T = 
4.2 K, in Faraday configuration, at i? = 14 T (:^ « 0.45), 
and for the circular polarization of light (correspond- 
ing to optical transitions of an electron in the excited 
spin state). Emission intensity / is plotted as a func- 
tion of both £in and Sont- In the inset we showed the 
following three cross-sections, (i) PL (emission) spec- 
trum is /(£out) for a fixed, high Sm. Although it is typ- 
ically measured for a much higher Sm, here we plot the 
data corresponding to resonant excitation of X^. (ii) 
PLE spectrum is /(fin) for a fixed, low fout, here corre- 
sponding to the recombination. For regular depen- 
dence of relaxation on PLE is an indirect measure 
of absorption, (iii) Diagonal cross-section /(fout = fin), 
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FIG. 10: (color online) Polarized photoluminescence (PL) 
spectra measured in a symmetric w = 15 ran GaAs quantum 
well with hole concentration g — 1.5 ■ lO^'^ cm~^, in magnetic 
field B = 14 T, at temperature T = 1.8 K. 



usually difficult to detect due to strong reflection of the 
incident light from the surface. By smoothing and ap- 
propriate orientation of the surface with respect to the 
incident/reflected direction we were able to minimize this 
effect and record meaningful diagonal spectra. 

In Fig. two strong peaks correspond to two optically 
active trions, and X^. Other, minor features are 
artifacts (note a sizable 0.12 meV step in £in). The sup- 
pressed intensity of the X (whose position relative to the 
X+'s is anticipated from the following Fig. I1U|I relative 
to the bright X+'s results from a rather large hole con- 
centration. On the other hand, the strongest X^ peak 
along the diagonal confirms earlier theoretical prediction 
for the trion oscillator strengths.'^" 

In Fig. 1101 we plot a pair of polarized cr^ and (t+ PL 
spectra /(£out)- They were recorded in the same sample 
and at the same magnetic field as PLE of Fig. ^ but 
for a higher excitation energy fi„ = 2.57 eV (wavelength 
488 nm) and at low temperature T = 1.8 K. Three trion 
states {X^, X^, and X^) along with the exciton are 
identified for the polarization. In the spectrum, 
the X^^ peak is weakened and X^^ disappears completely 
due to spin polarization. Additional low-energy peaks 
marked as AX^ = A'^ + X^ show recombination of a 
spin doublet {S — ^) ground state of a trion bound to an 
neutral acceptor located inside the well (cf. Fig.0J; such 
impurity-bound trions are discussed elsewhere^i 

The unambiguous assignment of the peaks was possi- 
ble from the analysis of the field evolution of the spectra, 
from B = beyond i? = 14 T, presented in a separate 



publication.^^ As noted by Glasberg et aZ.,A5i the Zeeman 
splitting of different X and X~^ states is very different 
(because of the wavevector dependence of the Lande g- 
factor for the holes). We marked these splittings with 
color horizontal bars in Fig. llOt b). The large value of 
1.6 meV for X^ (and AX^) compared to only 0.7 meV 
for X and 0.6 meV for X^^ is related to the occupa- 
tion of a strongly repulsive zero-angular-momentum hole 
pair state, characteristic of the singlet trion (and of dou- 
blet AX^). The analysis of the pair-correlation function 
shows that this pair state is not occupied also in 
(despite triplet spin configuration). This lets us expect 
similar X^ and X^ Zeeman splittings, even though X^^ 
is only detected in one polarization. 



Knowing the difference between X and X+ Zeeman 
splittings is necessary for a meaningful comparison of the 
trion binding energies A with the calculation of Fig.l^Jb) 
in which the Zeeman energy was ignored. The Coulomb 
binding energies A are extracted from the experimen- 
tal PL spectra by comparing the average X and X^ 
energies measured in both polarizations,^^ in Fig. llOf b) 
marked by dots on the Zeeman bars. In this way, we find 
A+ = A^^ « 1.4 meV and A^ w 0.6 meV. Compared 
to these values, the x = 1 calculation of Fig. [2Ib) pre- 
dicting A = 1.2, 1.35, and 0.4 meV, respectively, slightly 
underestimates the binding of all three states. The slight 
discrepancy could result from including only two low- 
est quantum well subbands and five lowest LLs, assum- 
ing equidistant LLs also for the holes, and ignoring the 
light- hole/heavy- hole mixing, all likely to enhance bind- 
ing. The numerical tests indicate that neither higher 
subbands (s > 2) nor the variation of LL spacing (be- 
yond n = 1) play a role, but the the n < A restriction 
indeed has a noticeable (< 0.2 meV) effect on A. 



The spectra in Fig. are analogous to Fig. ^| but 
they were recorded on an electron gas and involve nega- 
tive trions. The sample is also symmetric w — \b nm 
GaAs quantum well, but the concentration is slightly 
higher, g — 2 ■ 10"'^^ cm~^, so we had to use a stronger 
field, B — 22.5 T to keep a sufficiently low filling fac- 
tor, V « 0.4. In both polarizations, the detected tri- 
ons correspond to those of Fig. ^| The weak X peak 
marked for cr~ can be identified more convincingly when 
intensity is plotted in the logarithmic scale. The Zeeman 
splittings found in this system are 1.8 meV for X~ and 
X^^ (the same value was hence assumed for X~^^), and 
1.4 meV for X. Weaker variation of the splitting supports 
its attribution to the occupation of two-hole states. Us- 
ing these splittings, we find the following Coulomb trion 
binding energies: AjT « 2.1 meV, A^J^ 1.8 meV, and 
A^ « 0.65 meV. Again, calculation of Fig. Ufa) pre- 
dicting A = 2, 1.4, and 0.3 meV, respectively, slightly 
underestimates all these values. 
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FIG. 11: (color online) Polarized photoluminescence (PL) 
spectra measured in a symmetric w = 15 nm GaAs quan- 
tum well with electron concentration q = 2 ■ 10^^ cm~^, in 
magnetic field B = 22.5 T, at temperature T = 1.8 K. 

V. CONCLUSION 

We have carried out exact numerical diagonalization 
of realistic e + h and 2e + h hamiltonians (including spin, 
Coulomb interactions, and mixing of LLs and quantum 
well subbands) on a Haldane sphere. The parameters 
used for illustration are adequate for a symmetric 15 nm 
GaAs quantum well in a magnetic field i? = 10 T. 

Calculation of trion binding energies has been consid- 
erably advanced by taking the actual subband wavefunc- 
tions for the integration of Coulomb matrix elements and 



by including five LLs and two quantum well subbands for 
each electron and hole in the exact diagonalization. 

From the full energy spectra we have calculated the 
density of states (DOS). The main difference between 
DOS of e + h and 2e + h, representative for the exci- 
tons with and without the presence of additional free 
electrons, is the emergence of discrete bound trion states 
below the excitonic tails and redistribution of DOS away 
from the LL peaks and filling the gaps between them. 
The effect of an impurity on the e + h DOS is also stud- 
ied as a function of its charge and position in the well. 

For the full spectra of eigenstates, we have calculated 
the optical oscillator strength il. The optical density of 
states (ODOS) of 2e + h shows no obvious correlation 
with DOS. It is fairly complicated, with a great number 
of strong transitions labeled by LL indices of the recom- 
bining e-h pair and of the left-over electron. However, 
the e 2e + h absorption spectrum is far simpler. It is 
dominated by a series of LL peaks for the purely excitonic 
transitions vac —^e + h. These main peaks are neither 
shifted in energy, nor is their intensity noticeably sup- 
pressed or enhanced. The presence of (and interaction 
with) an additional electron shows in form of additional 
weaker peaks. Some of them are attributed to bound 
trion states (in the lowest and higher LLs). 

The numerical results have been successfully compared 
with PL/PLE experiments carried out on a 2D hole and 
electron gases. In particular, absorption and emission of 
the whole family of both negative and positive trions in 
the lowest LL has been observed. 
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